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ABSTRACT 

It is of great important to study pulsar beam shapes if we are concerned pith emission 
theories and pulsar birth rate. Both observations o and/or the ICS modelH2l show that 
different emission components are emitted from different heights. The relative longitude phase 
shifts due to different heights of emission components and to the toroidal velocity of electron 
are considered in this paper. Several possible observational features arising from the phase shift 
effects are presented. The emission beams may not have circular cross sections although the 
emission regions may be symmetric with respect to the magnetic axes. 

Subject headings: pulsars — polarization — radiation mechanisms 



1. INTRODUCTION 

Soon after the discovery of Pulsars, the rotating 
vector model (RVM) was proposed to well explain 
the linear polarization characteristics of pulsar ra- 
dio emission, in which Radhakrishnan and Cooke™ 
assumed that relativistic particles stream out along 
dipolar magnetic fields of pulsars and the polariza- 
tion position angle of emitted radiation from these 
particles is correlated with the magnetic curvature. 



Theoretical extensions of t 
literatures (e.g., Ferguson 



Le RVM have appeared in 
J, Blaskiewicz et. al.o, 



Hibschman Sz Aronsa), where some detailed physi- 
cal factors, such as the special relativistic effects, the 
polar-cap current flow, have been included. 

However the RVM fails in attempting to explain 
the position angle jumps observed in mean pulses as 
well as in individual pulsesO. This leads to the as- 
sumption that the radio beams of pulsars contain two 
nearly orthogonal polarization modes, the superposi- 
tion of which results in the non-S-shaped position an- 
gle variation in mean pulses and two ~ 90° separated 
distributions of position angles in individual pulses 
(see the paper by Mckinnon Sz StinebringB and ref- 
erences therein). Nevertheless, there is still a possi- 
bility that modified RVM, without inclusion of the or- 
thogonal modes, can reproduce the observed position 
angle "jumps" if we consider the linear depolariza- 
tion of pulsar beams and the observational unperlain- 
ties originated from observation polarimetersEI B . In 
view of the observational fact that radio emission of 
pulsars may have three components (core and two 



cones 



. t which are emitted from different heights above 
pulsar surface®, Xu et al.@ found that the posi- 
tion angles of mean pulses jump at certain longitudes 
where the linear polarization is approximately zero if 
the retardation effect, which causes the beam centers 
to be shifted between each other, can significantly 
conduce toward the linear depolarization. In this pa- 
per, we further study this issue arising from differ- 
ent emission heights which are calculatecLaccording 
to the inverse Compton scattering modelH in which 
core and conal components can be naturally repro- 
duced. 

Besides the polarization effect discussed above, dif- 
ferent emission heights may also have notable conse- 
quence to beam shapes. Many authors investigated 
observationally the pulsar beam shapes, the previ| 
ous results of which are summarized below. JonesEJ 
find that the latitudinal radius of beam cross section 
could be 2.5 | _limes the longitudinal one. Narayan & 
Vivekanandlimodify this ratio to be 3.0. Lyne & 
Manchester's" conclusion is that the shape of 
emission beam is approximately circular. Biggs'! 
proposes that therejnay be some meridional compres- 
Wu & ShenO find that, for pulsars with short 



sion. 



period, the radius of the emission beam at latitudinal 
direction is larger. However, the latitudinal radius is 
smaller for pulsars with longer period. Rankin, there- 
fore, change her morphological beam shape from an 
elliptical cross sectiono to a circular oneO . Our nu- 
merical results in section 3 may be checked by such 
observational study in the future. We find that the 



email: rxxu@bac.pku.edu.cn, xjw@bac.pku.edu.cn, gjn@pku.edu.cn 
2 This work is supported by NSFC (No. 19803001), by the Special Funds for Major State Basic Research Project, and by Doctoral 
Program Foundation of Institution of Higher Education in China. 



2 



emission beams of pulsars with small rotation periods 
do not have circular cross sections even if the emission 
regions are symmetric with respect to the magnetic 
axes if we take into account the toroidal velocity due 
to rotation. 

This paper is organized as below. After inspect- 
ing various aspects of the retardation effect caused 
by different emission heights and their possible ob- 
servational consequences in section 2, we calculate 
the emission beams from relativistic particles with 
poloidal (along magnetic field lines) and toroidal 
(perpendicular to field lines ) velocities in section 3. 
A conclusion, based on which we can know the phase 
shift between two emission units at arbitrary posi- 
tions, is proved geometrically in Appendix. The pa- 
per is summarized in section 4. 

2. RELATIVE LONGITUDE SHIFT DUE TO 
DIFFERENT HEIGHTS 

It is generally believed that pulsars' radio emission 
is originated from the region near the last open field 
lines. We assume in this section that the particles in 
the last-open-field-line region, which are responsible 
to observed radiation, have only poloidal velocities 
along the filed lines (i.e., the toroidal velocity is ne- 
glected), and are accelerated symmetrically respect 
to the magnetic axis. The emission region of each 
beam component is thus also symmetric respect to 
the axis. According to Appendix, the phase shifts of 
all emission elements of a certain beam component is 
the same since the values r • no of those emission ele- 
ments are the same, where r is the vector position of 
an element, no the emission direction. Therefore the 
emission beams should have circular cross sections in 
this case. 

2.1. Emission heights and beam shift 

We know from Appendix that an emission ele- 
ment at vector position r is equivalent to a virtual 
one whose longitude phase shifts to an earlier value 
8(p(r) = Qr ■ no/c. We now calculate 84>(r) func- 
tion of r, based on the dipole geometry of magnetic 
field configuration. 

For a dipole field line, we have 



A^sin 2 9, 



(1) 



where r = |r|, G is the polar angle, P the rotation 
period, c the speed of light, and A a parameter which 
characterizes the sorts of field lines (A ^ 1, A = 1 for 
the last open filed lines). Based on Eq. (jlj) , the angle 
between the direction of magnetic field at vector 
position r and the magnetic axis [i reads, 



cos ©, 



2\cP-6irr 



y/\cP(4:\cP-6Trr) ' 
2-3 sin 2 6 



(2) 



Therefore, assuming that elementary emission is 
along magnetic field lines and according to Eq.pH2[) 
and Appendix, we obtain the phase shift 8<j> as a func- 
tion of r 



4"7rr / XcP — 2nr 



XcP V 4AcP - 6vrr ' 



(3) 



and the phase difference A<f> between two emission 
beam components with heights r a and r^ is 



= <ty(r a ) - 8<j>(r h ). 



(4) 



We note, from Eq.Q, that 8(f) is only a func- 
tion of rotation period P and distance r, but 
is not relevant to magnetic inclination angle a 
and/or impact angle (3. The calculations of phase 
shift 8(f) for pulsars with periods P = 0.03s, 
0.05s, 0.1s, 0.5s, and Is are shown in Fig.||. 
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Fig. 1. — The phase shifts 8<f> for pulsars with P = 0.03s, 
0.05s, 0.1s, 0.5s, Is, based on Eq.([]). The emission height 
(spherical coordinate r) is from 30km to 10 3 km. A = 1. 

We see from Fig.|l| that the retardation effect arising 
from different emission heights is less important for 
pulsars with longer periods P, but can not be negli- 
gible when P is small. For instance, the phase differ- 
ence between components with heights of 30km and 
300km could be over 10° if P = 30ms. But the phase 
difference can only be about 1° for emission compo- 
nents with heights 50km and 500km if P = 0.5s. 



2.2. The relationship between observed frequency 
and shifted phase 

Theoretically, we can know the emission heights of 
beam components at a certain observation frequency 
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in the inverse Compton scattering modelMiJ , based 
on which the three components (core, inner cone, 
and outer cone) can be understood naturally. The 
frequency v of radio wave emitted at position r (or 
height r — R, R is the pulsar radius) can be got- 
ten by following Eq.([|-[7]). Low frequency electro- 
magnetic waves with frequency are supposed to be 
produced near_pulsar surface due to RS-type vacuum 
gap sparkingH. These waves are assumed to prop- 
agate nearly freely in highly inhomogeneous plasma 
in pulsar magnetospheres and inverse Compton scat- 
tered by the secondary particles with Lorentz factor 
7j_turning into radio waves observed with frequency 



1.5YW1 



1-7- 



' cos 1 



(5) 



where the incident angle Oi is the angle between the 
wave vector of low frequency wave and the direction 
of electron moving (along magnetic fields in the ap- 
proximation of this paper), which can be calculated 
to be 



cos 6>i 



2cosG + (R/r)(l -3 cos 2 9) 



y/(l + 3 cos 2 9)[1 - 2{R/r) cos 6 + {R/r) 2 \ ' 

The Lorentz factor of relativistic electron should de- 
crease with height 3 , which is suggested in the follow- 
ing form 113 



7 = 7 exp[-£- 



R, 



R 



(7) 



where 70 is the initial Lorentz factor near pulsar sur- 
face, £ a parameter reflecting the extend of the energy 
loss. 
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Fig. 2. — The phase differences between beam components, 
as functions of observation frequency, in the ICS model. The 
upper two lines (solid and short dashed) are the phase differ- 
ences of outer cone and inner cone, while others are of inner 
cone and core. The pulsar periods are chosen to be 0.1s and 
0.05s, and other parameters to be 70 — 3 x 10 3 , £ = 3 x 10 -2 , 
R = 10 4 m, and v = 10 5 H Z . A = 1. 

The relation between observation frequency v and 
position r can be found from Eq.(jl], ^-0). For a cer- 
tain v we may find three solutions for r: ri , r<z, T3, be- 
ing the positions of core, inner cone, and outer cone, 
respectively. The phase differences of beam compo- 
nents of outer and inner cones, and of inner cone and 
core, are calculated for observation frequency up-to a 
few GHz in the ICS model, which are shown in Fig.||. 

We see from Fig.|2| that, on one hand, at lower fre- 
quency v, the phase differences Afo = 8<j) (r^)— Sfifa) 
of outer and inner cones are notable. For example, 
A<p2 = 30° when v = 2 x 10 7 Hz. A</>2 decreases 
monotonously with the increase of observation fre- 
quency v. On the other hand, the phase differences 
A0i = 5(j){r2) — 6<f)(ri) of core and inner cone compo- 
nents are small at lower v. A<^>1 = 0.05° contrary to 
A</>2) A</>i increases monotonously with the increase 
of v. The observational tests and/or consequences of 
these features need further investigations in the fu- 
ture. 

2.3. The possible distribution of singular points on 
the celestial sphere 

There may be some points on the observational ce- 
lestial sphere where the linear polarization intensity 
L = 0. We call such points as "singular points". 
The polarization position angle should jump 4 exactly 
90° when a fine of sight goes across one of the singu- 
lar pointsUJ 0. Therefore the distribution of singular 
points on the celestial sphere would determine the 
variation properties of position angle. 

Generally, it is hard to calculate the singular point 
distribution since there are many factors, heing not 
well understood, to cause depolarization!-' ___3 of pul- 
sar beamed radiation. However, if depolarization is 
mainly due to the incoherent superposition of emis- 
sion components at different heights, and if there are 
only two components, we can get the singular point 
distribution easily in this special case. 

In this case the sufficient and necessary conditions 
that should be satisfied at a singular point are 



tp2 



±vr/2 
L 2 



(8) 



where V'i and L\ (i = 1, 2 denotes the first and second 
components, respectively) are the position angles and 



Electrons will lost their kinetic energy, via, e.g., the curvature radiation and/or the inverse Compton scattering of these electrons 
off thermal X-ray photons from pulsar surface, when they are moving along magnetic field lines. 
4 The position angle may jump approximately 90° if a line of sight goes nearby a singular point. 
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linear polarization intensities of i components, re- 
spectively. The first condition, i.e., tpi — ip2 = ±tt/2, 
is thus only a necessary one for a singular point. Con- 
sidering the identity tan(^-^ 2 ) = g^gfe , we 
can re-write the necessary condition below 



tan tpi tan ip2 



0) 



The possible distribution of singular points on the 
celestial sphere is calculated for this special case 
based on Eq.Q, and is shown in Fig.||. We assume 
the phase difference between the centers of the com- 
ponents due to their different emission heights to be 
A(p = 15°, and the inclination angle a = 60° or 30° in 
our calculation. We find from Fig.|3] that the distribu- 
tion become "flatter" when the inclination angle a is 
smaller, i.e., larger a may be more favorable for a line 
of sight to go across a singular point, thus increasing 
the possibility of position angle curve to "jump". 
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Fig. 3. — A possible distribution of singular points on the 
celestial sphere assuming that there are only two emission com- 
ponents, and that the linear depolarization arises only from the 
incoherent superposition of those two components. The phase 
difference between the components is assumed to be 15°, and 
the inclination angles are chosen to be 30° and 60°, respec- 
tively, in the computation. 




Fig. 4. — A sketch of two coordinate systems, we assume 
the /^-coordinate, in which the magnetic axis fj, is chosen as 
the z-axis, turns clockwise an angle a (the inclination angle) 
along the x-axis of f2-coordinate. A position vector r can be 
represented by {r, 8, cf> } in Sl-coordinate, or by {r, ©,$ } in 
/x-coordinate. rie is the vector unit of magnetic field. 

In the above discussion, only the electron veloc- 
ity along magnetic field line (the poloidal velocity) 
is considered. However, pulsar rotates, an elec- 
tron does have a toroidal velocity fi x r (see FigQ). 
Actually, the velocity of an electron in the observer 
rest frame can be obtained by the relativistic veloc- 
ity transformation of the poloidal velocity in the ro- 
tation frame. The inclusion of such toroidal velocity 
has at least one implication: an emission unit would 
shift forward a phase with respect to the magnetic 
direction riB since a relativistic electron emit pho- 
ton at the direction of its moving. In this section, 
we calculate the beam shapes modified by inclusion 
the toroidal velocity, basexLon the derived emission 
heights in the ICS model t3. In the following com- 
putation, two phase shift effects arising from differ- 
ent emission heights and from toroidal velocity are 
included. According to Fig.[|, one can find the rela- 
tionship between {#,</>} and {0,<1?}, 




cos a cos — sin a sin sin $ 

cos a sin <E> sin 0+sin a cos B 
sin 9 

cos <3> sin © 



(10) 



which is useful below. 

The linear rotation velocity at a position vector r 
is v. 



3. EMISSION BEAMS IN THE ICS MODEL 



v = \ft X r| 

= Qr yjl — (cos a cos — sin a sin sin $) 2 

(11) 

Thus, by the relativistic velocity transformation, one 
finds two components, u\ (in the flxr direction) and 
U2 (in the [(fi x r) x ne] x (ft x r) direction), of the 
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velocity of an electron in the observer rest frame, 



u 2 



'• sin a sin © cos $+i>V 1+3 cos 2 



\/l+3cos 2 0— vy/l— 7 — 2 sin a sin 9 cos $/c ' 
c^/ (1— 7~ 2 )(1— u 2 /c 2 )(l+3 cos 2 0— sin 2 a sin 2 ©cos 2 $) 
Vl+3cos 2 0— v^l— 7 — 2 sin a sin ©cos 3>/c 

where 7 is the Lorentz factor of electron in the ro- 
tation frame. Therefore the angle 5 between riB and 
the emitted photon is, 
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-If 



(13) 

which is also the phase shift due to toroidal velocity. 

Combining the phase shifts of 5(f) (Eq.@) and 5tp 
(Eq.([l3|)), we can calculate the emission beams in the 
ICS model. We use the results in Fig.|| to decide the 
emission heights. The radio emission is assumed to 
originated from a region between magnetic field lines 
denoted by A = 1 and that by A = 1.8. The radius 
of the stellar surface of the feet of the field lines cor- 
responding to A = 1.8 is ~ ~ O.750 p , where 

P = sin -1 \J^rjT is the radius of polar cap. The cases 
of rotation period P = 0.1s and 0.05s, and observa- 
tion frequency v = 10 9 Hz are presented in FigjB] and 
Fig|| The parameters in the simulation are listed in 
Table 1. 

From Fig|3] and Fig.[|, we see the emission beams 
are clearly not symmetric although the emission re- 
gions of these components may be symmetric respect 
to the magnetic axes. The width of the first pulse 
in a conal component would be much smaller than 
that of the second one. Besides the clear phase shifts 
between the components, the emission beams may 
not have circular cross sections when the rotation pe- 
riod is small (see Fig.^). The pulses of impact angle 
j3 > arrive earlier than that of /3 < 0. This is un- 
derstandable. Since the emission region of (3 > is 
farther from the rotation axis than that of /3 < 0, 
the phase shift effect thus is more remarkable due to 
larger toroidal velocity. 



Fig. 5. — The emission beams of core, inner cone, and outer 
cone simulated in the ICS model, with the inclusion of the 
phase shift effects both due to different emission heights and 
to toroidal velocity. The inclination angle a — 30°, the rota- 
tionperiod P — 0.1s. Other parameters are the same as in 
Fig! 
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- Same as in FigJB|, but for P = 0.05s. 
4. CONCLUSION AND DISCUSSION 



A detail consideration of the relative longitude 
phase shift due to different emission heights of three 
components is presented in this paper. Our main 
conclusion can be summarized below. 

1, An emission unit along direction ni(0o,0i) at 
an arbitrary position ri arrives earlier than another 
emission unit along direction 112(80,4*2) at another 
position r2, the phase difference between those two 
elements is <f>\ — 4>2 + f2(ni ■ ri — 112 • ^2)/c (if this 
value is positive), where O is the rotational angular 
velocity. 
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Table 1 

The parameters* chosen and calculated in the ICS model for the computations of 

EMISSION BEAM COMPONENTS 



*P - rotation period; v - observation frequency; r; - emission height of component i; i= 1, 2, 3 denotes core, 
inner cone, and outer cone, respectively; 9^ - angular radius of component i; A(fri and A<fo are the phase 
differences of core and inner cone, and of inner and outer cones, respectively. 



P(s) 


zv(Hz) 


ri(m) 


r 2 (m) 


r 3 (m) 


0^1 


9^2 


9/13 


A0! 




0.1 


10 9 


18820 


81900 


483820 


5°.4 


11°.3 


28°. 1 


0°.76 


4°. 74 


0.05 


10 9 


20695 


56040 


697000 


8°.0 


13°. 2 


50°. 5 


0°.85 


14°. 71 



2, The phase shift effect is less important when the 
rotation period P > 0.5s. 

3, The phase shift values, as a function of observa- 
tion frequency v , are presented in the ICS model. As 
v increases, the shift between outer and inner cones 
decreases, while the shift between inner cone and core 
increases. 

4, Larger inclination angles may be favorable for 
the appearance of polarization position angle jumps. 

6, The emission beams are not symmetric in the 
ICS model, which is achieved by the inclusion of 
the phase shift effects both due to different emission 
heights and to toroidal velocity (fi x r). 

In the simulation of the beam shapes of this paper, 
we just take the symmetric emission region in the 
ICS model for simplicity. However, the emission re- 
gion should actually not be symmetric respect to the 
magnetic axes when we consider that, for a certain in- 
verse Compton scattering process, the low-frequency 
wave is emitted at an earlier sparking point, and that 
electrons are not moving exactly along magnetic field 
lines. How about the emission beams and polariza- 
tion properties when these issues are included? A 
further study is needed in the future. 



acknowledgements: We thank our pulsar group for 
discussions. 



that: first we find a virtual equivalent 5 at origin for an 
emission unit at arbitrary position; then we compare 
the phases of the virtual equivalents. The phase dif- 
ference between two virtual equivalents is actually the 
phase difference observed of those two corresponding 
emission units. Our conclusion can be expressed as 
following proposition: 

An emission unit along a direction no(#0;^o) a t an 
arbitrary position r(9, <f>) is equivalent observationally 
to a virtual emission unit along direction nQ(#o,0o) 
at origin, where (fr' = (fro + Sir, r = r • uq/c, f2 and 
c are the rotational angular velocity and the speed of 
light, respectively. 



n' (e o ,cfO 




Appendix 

Observation shows that emission components (core, 
inner cone, and outer cone) could come from different 
heightso. The longitude phase shifts between those 
components are thus inevitable!! . We now try to 
determine this shift value between two emission com- 
ponents. Our measures to deal with this question is 

5 We call two emission units are equivalent only if an observer at infinity can not distinguish them, i.e., the observer detects same 
intensity at same time for those two emission units. 



Fig. 7. — A sketch geometrical figure for proving the conclu- 
sion to find an equivalent emission unit at origin. An emission 
along no at point E (position vector r) is equivalent to the 
virtual one along n' at origin. 

We are to prove this below. As in Fig.[/], we erect co- 
ordinates o-xyz in the observation rest frame, where 
z-axis is along rotational axis (O), the origin is at 



the center of pulsar. We are considering the emis- 
sion along no at an arbitrary point E. A plane ABC 
is created so that no is perpendicular to this plane. 
Because all emission units along the direction no at 
those points which are on the ABC plane should be 
observed simultaneously, the emission along no at 
point 6 F is equivalent to the emission at point E if 
the intensities are the same. We assume the time 
t = when a photon is emitted at point E along no- 
A photon emitted at point O along no when t = — r 
(r = |OF|/c) arrives at point F when t = (still 
along no). However, at t = 0, the emission at point 
O is along n' (6>o, (fto + Qr) because of rotation. There- 
fore, an emission along no at point E with position 
vector r is equivalent to the virtual one along n' at 
origin O. 

Now we calculate the value r. In the right-angled 
triangle OEF, |OF| = rn . We thus have r = t-toq/c. 
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Point F is the intersection of ABC plane and its perpendicular line through point O. 



